High-order Ho multipoles in H0B2C2 observed with 
soft resonant x-ray diffraction. 
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Abstract. Soft resonant x-ray Bragg diffraction (SRXD) at the Ho M4.5 
edges has been used to study Ho 4/ multipoles in the combined magnetic 
and orbitally ordered phase of HoB2C2- A full description of the energy 
dependence for both cr and tt incident x-rays at two different azimuthal angles, 
as well as the ratio la- /In as a function of azimuthal angle for a selection 
of energies, allows a determination of the fiigher order nmltipole moments of 
rank 1 (dipolc) to 6 (liexacontatetrapole). The Ho 4/ multipole moments 
have been estimated, indicating a dominant hexadecapole (rank 4) order with 
an almost negligible influence from either the dipole or the octupole magnetic 
terms. The analysis incorporates both the intra-atomic magnetic and quadrupolar 
interactions between the 3d core and 4/ valence shells as well as the interference 
of contributions to the scattering that behave difTerently under time reversal. 
Comparison of SRXD, neutron diffraction and non resonant x-ray diffraction 
shows that the magnetic and quadrupolar order parameter are distinct. The 
(00 1 ) component of the magnetic order exhibits a Brillouin type increase below 
the orbital ordering temperature Tq , while the quadrupolar order increases more 
sharply. We conclude the quadrupolar interaction is strong, but quadrupolar 
order only occurs when the magnetic order gives rise to a quasi doublet ground 
state, which results in a lock-in of the orbitals at Tq. 



1. Introduction 

Rare earth based materials are well known for their applications as strong permanent 
magnets such as the intermetallics Nd2Fei4B and SmCo5. The large orbital 
momentum of the rare earth 4/ electrons is essential for the strong magnetic anisotropy 
of these hard magnets. The 4/ valence electrons are localized and they may retain 
spin and orbital degrees of freedom down to very low temperature, which can also 
result in novel ordered states such as the Kondo effect, heavy fermions, complex (sine. 



High- order Ho muUipoles in H0B2C2 



2 



spiral, triple-k) magnetism, and orbital order [HIS]- 4/ orbitals are described to first 
order by their quadrupole moment although higher order moments may also order, 
and this can occur independently from the magnetic order. 

Many functional properties such as multiferroicity, magnetoresistivity and 
superconductivity arise from competition between distinct order parameters. The 
frustration of normally dominant interactions can make higher order terms relevant. 
For example, ordering of the electromagnetic anapole has been observed recently 
in CuO [3] and may be relevant to both its multiferroic properties and the 
superconductivity of related compounds. Ordered states of higher rank multipole 
moments appear in selected compounds, such as octupoles in Cei-j^Laa^Bg [H [S], 
hexadecapoles (or triakontadipoles) in Np02 [3 |7], and hexadecapoles in URu2Si2 
[3 Hj. We have investigated the interplay between the magnetic and quadrupolar 
order of the 4/ electrons in H0B2C2 using resonant x-ray diffraction (RXD). 

RXD, with its element specificity and enhancement of the scattered intensity by 
up to three orders of magnitude, has been shown to be an invaluable tool in the 
investigation of electronic properties of materials, especially in cases where there are 
multiple contributions to the scattered intensity such as spin, charge, and orbital 
orders [TD]. In materials where the electrons are largely localized, it is advantageous 
to employ a local approximation such as an expansion in spherical multipole moments 
Y^, and their spherical tensor representation (where x < L and —x < q < x). Such 
spherical multipoles are parity even and have the definite time signature (— 1)"^. We 
refer to time-odd multipoles as magnetic and time-even multipoles as orbital. X-ray 
diffraction has a convenient description in terms of spherical tensors, which includes 
information on the scattering geometry and polarization, and the scattering factor 
becomes a simple tensor product sensitive to properties of the x-ray beam and the 
different spherical components of the electron density at the rare earth site. 

RXD is particularly successful in the investigations of exotic orders, and has for 
example demonstrated the existence of ordered multipole moments of up to rank 6 
(the maximum allowable for a rare earth atom) in DyB2C2 [HI [H]. While RXD 
studies of rare earth elements typically probe the 2p-5d dipole (El) transition and the 
often much weaker 2p-4f quadrupole (E2) transition, soft resonant x-ray diffraction 
(SRXD) probes the 3d-4f dipole transition, and the strong interaction between the 4/ 
valence and 3d core shells can give rise to complex spectra. 

At room temperature, H0B2C2 crystallizes in the tetragonal space group PA/mbm 
(see Figure [1]) and undergoes an antiferromagnetic (AFM) phase transition at = 
5.6 K to the so-called Phase IV, and a further transition to a combined AFM and 
antiferroquadrupolar (AFQ) ordered phase at Tq = 5 K (termed phase III) which 
has been confirmed directly with RXD at the Ho L edges [13] and also indirectly 
by ultrasonic measurement of the elastic constants [14]. Below Tn the magnetic 
structure is incommensurate with ordering vector [1 -I- 6,6,6'] with 6 = 0.11 and 
6' = 0.04. The magnetic structure undergoes a reorientation at Tq to an AFM order 
having propagation vectors fci = [100] and k2 — [01^] together with a ferromagnetic 
component ^3 = [000] and a weak component with k4 — [00|] 15 . 

This alternating arrangement of tilted moments is postulated to arise from 
competition between the magnetic exchange and coulomb type orbital interaction 
|16| . It is remarkable that the quadrupole transition occurs at a temperature lower 
than Tn , as typically the onset of magnetic order removes the degeneracy of any low- 
lying crystal field levels. Quadrupolar transitions are most often associated with the 
splitting of a quartet (or pseudo-quartet) into two doublets, which is ruled out if all 
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Figure 1. Magnetic and crystalline structure of H0B2C2 according to Ref. 
|16l . Arrows indicate the direction of the Ho 4/ magnetic moments in the low 
temperature Phase III. 



low-lying levels are separated into singlets by magnetic order. The reverse situation, 
where the doublets formed by quadrupolar ordering are again split into singlets by 
magnetic order, has been observed in several compounds including DyB2C2. Based on 
the release of entropy equal to Rln3 at the successive transitions it has been suggested 
[17j that the Ho ground state is a quasi-triplet consisting of a singlet and a low lying 
doublet, and that the doublet splits at Tn. In this scenario, below Tn one of the levels 
of the doublet approaches the singlet forming a quasi-doublet which then gives rise to 
the quadrupolar transition at Tq. Although many investigations of this material cite 
an influence from quadrupole order or quadrupolar fluctuations |18j . measurement of 
lattice constants has inferred the presence of an octupolar order parameter |14j , which 
together with the observation of high order multipoles in DyB2C2 suggests that such 
quantities could be observed in the low temperature Phase III. 

In the present work, we performed an SRXD study of H0B2C2 in phase III 
(combined AFM + AFQ) and collected neutron and non-resonant x-ray temperature 
dependences for reflections relevant to the orbital ordering. The core hole splitting 
formahsm previously applied to the AFQ and AFM -I- AFQ phases of DyB2C2 pT|[T2] 
is employed to interpret the energy and azimuthal dependence. In sections 2 and 3 
we give an overview of the experimental details and results, followed in section 4 by 
a theoretical description of the resonant scattering amplitude in terms of the state 
multipoles of the 4f shell. In Section 5 we discuss the results for the magnitude and 
sign of the quadrupolar and magnetic intra-atomic interactions, and the estimations 
of the strength of the high order multipole moments. We also discuss anomalies in the 
temperature dependence that suggest the magnetic and quadrupolar order is distinct 
below Tq. Section 6 concludes our paper with a summary of findings. 

2. Experimental Details 

A H0B2C2 single crystal was grown by the flux-pulling method using an arc-furnace 
with four electrodes and cut and polished with (001) perpendicular to the sample 
surface. The orbital ordering (00|) reflection was recorded with tt and cr incident 
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Figure 2. Scattering geometry used in tiie SRXD experiment indicating the 
polarization of the incident radiation, sample orientation and momentum transfer 
q = k - k'. 

polarization at the Ho M4^5 edges of H0B2C2 at the UE46-PGM-1 beam-hne at 
BESSY II synchrotron in Berhn, Germany [TO]. Figure [2] shows the geometry of 
the experiment, and a zero degree azimuthal angle corresponds to alignment of the 
crystallographic b-axis in the scattering plane. The sample was cooled to T = 3.5 K in 
a liquid helium flow cryostat and a thin aluminum thermal shield was used to maintain 
the temperature. Beam heating of the sample was observed which was minimized by 
reducing the beam intensity to a point where additional reduction did not change 
the measured ordering temperature. The energy dependence of the absorption was 
determined from the FWHM of 6* — 26* scans performed at each energy. The intensity 
of the (00 1) reflection was recorded as a function of both incident polarization 
and energy. At selected energies the temperature and azimuthal dependence were 
measured. The sample size was larger than the dimensions of the beam, making it 
difficult to compare absolute intensities as a function of azimuth. Consequently it is 
the ratio la / In that is analyzed in this work. 

In addition, the temperature dependence of the (00 1) reflection was measured 
with neutron diffraction (ND) on the WOMBAT high intensity neutron powder 
diffractometer at the OPAL reactor in Sydney, Australia using 2.995 A neutrons 
[50]. Further, the temperature dependence of the (01^) reflection was measured 
with non-resonant x-ray diffraction (XRD) at SPring-8 on the BL19LXU beamline 
[21] . Measurements were taken using 30 KeV x-rays and a Ge solid-state detector to 
eliminate higher-order harmonics. 

Heat capacity was measured using a Quantum Design MPMS at Toyama 
University, the results of which are presented in Figure [3j showing the onset of Phase 
HI at Tq=4.8 K (sharp peak) and the onset of Phase IV at Tjv=5.5 K (broad peak). 
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Figure 3. Heat capacity measurements on H0B2C2, sliowing the onset of Phase 
III at Tq = 4.8 K (sharp peak) and the onset of Phase IV at Tjv = 5.5 K (broad 
peak) 

3. Experimental Results 

Figure |4] sliows tiie recorded absorption (a,b) and integrated intensity (c,d) of the 
space-group forbidden (00|) reflection of H0B2C2 in phase III at an azimuthal angles 
of (f —45° and ip =90° as recorded with SXRD. Also shown is the absorption corrected 
intensity (e,f) obtained by taking the product of the integrated intensity and the 
absorption. We draw particular attention to the qualitative difference between the a 
and TT spectra. The intensity of the central features labeled B and C at the M5 edge 
observed for ip — 45° with cr incident radiation are much larger when compared to 
the other spectra, while the features A and D are largest for (p — 90° with tt incident 
radiation. This is indicative of the coexistence of two order parameters, namely the 
magnetic and orbital order, giving rise to interference. The significant extra intensity 
at (/? = 45° can be attributed directly to the orbital component of the scattering 
as will be demonstrated later. See also Appendix C. 

Figure [5] shows the intensity ratio 1^1 1-n as a function of the azimuthal angle 
ip compared to a single order parameter possessing 2-fold rotational symmetry (i.e. 
the 2/m Ho site symmetry). Such a cos^ 2(p dependence has been observed for the 
quadrupole E2 resonance {2p to 4/ transition) at the L edge [T3]. The azimuthal 
dependence at E= 1357 eV (corresponding to D in Figure 3) departs from this 
simple relationship, further confirming the interference of magnetic and orbital order 
parameters. 

Figure |6] compares the temperature dependence of SRXD recorded at E w 1347 
eV and E ss 1357 eV (labeled features B and D respectively in Figure [4]) for tt and 
cr incident radiation and = 45 ° and 90 ° with the recorded XRD intensity of 
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Figure 4. (a) and (b), absorption spectra of tiie (00 g) reflection in H0B2C2 
determined from 6 — 20 scans with incident a (vertical) and vr (horizontal) 
polarization at 3.5 K with an azimuthal angle of 45° and 90° respectively, (c) 
and (d), integrated intensity of same, (e) and (f) are intensities of (c) and (d) 
corrected for absorption. The marked positions A-D indicate the energies at which 
temperature and azimuthal dependence measurements were made. Note: All data 
are normalized to the central peak at the M5 edge in the cr channel, and the ratio 
between a and it channels is explicitly maintained. The intensity at the M4 edge 
is multiplied with a factor of 16 for clarity. 

the (0 1 11/2) reflection and ND intensity of the (00^) reflection. The temperature 
dependences measured with the three techniques are markedly different; while SRXD 
has a finite intensity above Tq both (00^) ND and (01^) XRD have zero intensity 
in this region. The temperature dependence measured with SRXD depends on 
polarization, azimuthal angle, and energy, and falls between the XRD and the ND 
temperature dependences. This will be discussed in sect, [sj 

4. Theoretical Pramework 

It has been demonstrated [11] [T^] that the unique shape of the energy dependence of 
the (00 i) refiection in DyB2C2 can be attributed to the splitting of the 3d core states. 
This splitting results in multiple interfering resonators and consequently adds structure 
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Figure 5. Azimuthal angle dependence of the intensity ratio 1^ / It: at selected 
incident x-ray energies corresponding to features (A through D, respectively) in 
energy dependence (Figure|4|. Solid lines are proportional to cos^ 2ip + constant, 
the expected azimuthal dependence based on the site symmetry. 
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Figure 6. Measured temperature dependence of the (00|) reflection recorded 
with SRXD at selected energies corresponding to specific features B and D (Figure 
[4J, for both (T and tt incident polarization, normalized at T = 3.5 K. The data is 
compared to the intensity recorded with XRD and ND. 
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to the energy profile. The sphtting of the core states is caused by a combination of the 
intra-atomic 4/ — 3d quadrupole interaction which splits the different \M\ levels, and 
the intra-atomic magnetic interaction which lifts the remaining degeneracy in ±|M|, 
where the bar refers to the quantum numbers that label the core-hole. The relative 
oscillator amplitudes are determined by the 4f wavefunction, which is characterized in 
terms of state multipoles (T^). 

The X-ray scattering amplitude, the square of which is proportional to the X-ray 
intensity, is formed from the product of two tensors |10j 

/ = E(-l)^x5,(2i. + 1) J2 DU-^ 7) E E ^-^'t:U^l)\^T^ 

with —K < Q < K . is constructed from the polarization vectors of the incident 

and scattered X-ray beams. The Wigner D tensor iDg^(Q!, /?, 7) rotates (which 
refers to local coordinate axes of the sample) , onto the coordinates of the experimental 
reference frame used for X^q with Euler angles a,/?, and 7. Our particular choice 
of principle axes gives D^^(a,^,7) = -Dq^(0, f ,7r - v?) = e'l'^D^^i^) The Ho 

site symmetry restricts q to ±1,±2. K = 2 corresponds to orbital order, if = 1 to 
magnetic order, and K = to charge order, of which the latter is absent at the (00 1) 
Bragg reflection. 

The sum over the describes the series of oscillators created by the splitting of 
the core state [TTl [T^ . The amplitude A^{J,M) of each oscillator is labeled by the 
total angular momentum of the core hole, J = |_for the M4 and ^ = | for the M5 edge, 
and by the magnetic quantum number —J<M<J. Regarding the denominator, E 
is the photon energy, Aj is the difference in energy between the degenerate 3dj shell 
and the empty 4f states and hT'j ^ is the lifetime of the intermediate state. e''( J, M) 
is the energy shift of the core levels due to the intra-atomic quadrupole interaction 
Qj and the magnitude of the intra-atomic magnetic interaction due to the unpaired 
4/ electrons i?int- 



e«( J, M) = [3A/2 _ J( J + l)]Qj + gjM . Hi^t (2) 

There are six resonant oscillators at the edge (M = ±|, ±|, ±i) and four at 
the M4 edge (M = ±|,±i). The oscillator amplitudes interfere and the branching 
ratio between the two edges is given by the purely real mixing parameter rj. A^ (J, M) 
is constructed from the structure factor of the chemical unit cell ^'^ and these two 
quantities are discussed in Appendicies A and B respectively. The Ho site symmetry 
dictates that A2 = — A?.2, and A^ = A^^, which follows from the properties of 
A} (J, M) and A\{J , M) are proportional to the dipole moment {Tl) , octupole moment 
(rf), and the triakontadipole moment (Tf) while A2{J,M) is proportional to the 
quadrupole moment the hexadecapole moment (Tj ), and the hexacontatetrapole 
moment (T^). The contriutions to (1) as a function of azimuthal angle and incident 
polarization are discussed in Appendix C. For further theoretical details, see [l2] . 

5. Discussion 

The theory detailed in section |4] is used to fit simultaneously the measured energy 
dependence for a and tt incident polarization and the measured azimuthal angle 
dependences via global optimization using a simple genetic algorithm method and 




Figure 7. Calculated energy dependence compared to normalized SRXD 
intensity at azimuthal angles of 90° (a and b) and 45° (c and d) for the AFM+AFQ 
phase at 3.5K. Presented are spectra measured with both tt (a and c) and cr (b 
and d) incident radiation, shifted upwards for clarity 
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Figure 8. Fitted azimuthal dependence of the intensity ratio Ict/Itt at energies 
corresponding to features A, B, C and D in Figure |4] shifted upwards for clarity. 
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a pattern search optimization of the genetic algorithm result. Although such a 
method cannot be guaranteed to locate the global minimum, by repeating the fitting 
process until the best solution has been obtained multiple times we obtain a statistical 
uncertainty {2a) of better than 10% for all parameters except the linewidth |3| 
which has an uncertainty of 20%. The best fit is shown in Figure [7] and Fi gure M for 
the energy and azimuthal angle dependence, respectively. The absorption correction 
is added to the calculated spectra so it can be compared with the as-measured data, 
and the oscillator positions are indicated in Figure [7] with arrows. The parameters 
obtained in the fit, and the associated angular distribution of the 4/ electronic charge 
density, are shown in FigurejOj The branching ratio r3/2 = 0.12 calculated for the free 
ion gives an excellent fit. We find that the intra-atomic quadrupole interaction 
Q5/2 = —0.09 eV, and the intra-atomic magnetic interaction |iJint| = 0.39 eV. Both of 
these values are largely independent of the choice of wavefunctions and linewidths, and 
the order of magnitude of Hint is in line with that expected from atomic calculations of 
3d-4/ exchange integrals [23]. It is important to note that whilst oscilators with ±M 
will be shifted in opposite directions by -ffint, there is no way to distinguish between 
them and so only the magnitude of i?int can be used as a fitting parameter, its sign 
being arbitrary. 

Variations in (T^) can be partially compensated by the line widths to yield a 
similar result. The linewidth for the oscillators Fj |^-^| vary from « 0.5 eV to 12 eV and 
vary substantially between time even and time odd contributions. The hexadecapole 
^2 is the dominant contribution and is normalized to 1, giving 5'} — 0.03, ^'2 — 0.36, 

= 0.07, and *f = -0.12, = 0.42. It is notable that the lower rank 
magnetic contributions (the octupole and dipole) are small relative to the time-even 
contributions. Whilst the actual ordered moment is given by (T^), the magnetic 
contributions (odd q) to the scattering are sensitive to {T^)' — (T^)" (where (. . .)' 
and (. . .)" represent respectively the real and imaginary parts of the atomic tensors, 
see Appendix A) which can potentially be small when the ordered moment is large, 
in the case where the real and imaginary parts are of similar size. This coincidence 
cannot be ruled out, as neutron diffraction studies have shown that the magnetic 
moment is large, although a substantial portion of the ordered moment propagates 
along directions other than the (00^) [55]. Additionally, the moment (J) measured 
with neutron diffraction is not precisely the same quantity as, for example, the (T^) 
magnetic dipole tensor which is a linear combination of (S), (L), and (3R(R • S) — S) 
[To] . Barring this accidental case of near-cancellation, the relatively small contribution 
of the ordered magnetic moments (of lower rank) can be attributed to the frustrated 
interactions of the system, as the magnetic and orbital exchanges strongly compete 
with one-another. 

For comparative purposes, we also plot a best fit in Figure [7] without the higher 
order multipoles (i.e. rank 5 and 6) that are not directly observable via any other 
means, and retain only terms of rank < 4 that are observable via an E2 resonance. 
The inclusion of higher rank terms significantly improves the description of the data 
in the pre-edge region, particularly where the intensity dips towards zero. 

The values of Q5/2 and |i?int| obtained by this analysis are smaller in magnitude 
and of the same sign as those obtained from the analysis of the AFM-I-AFQ phase 
of DyB2C2 where they were found to be —0.42 and 0.98 eV respectively [12 . The 
interaction between the 3d and 4/ shells is expected to be similar for both Ho'^+ and 
Dy'^^ as they contain 9 and 10 4/ electrons respectively, and the radial extent of the 4/ 
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shells are almost identical. This suggests both the magnetic and quadrupolar moments 
of the Ho 4/ shell are smaller than for the Dy 4/ shell in i?B2C2. The two materials 
have a similar magnetic structure with four magnetic sublattices characterized by the 
propagation vectors [10 0], [0 [10 ^] and [0 0] , and are inferred to have a 
similar pattern of orbital ordering on this basis, although H0B2C2 has a much lower 
ordering temperature |161 128j . The magnetic structures are compared in I28[ and the 
component of k along (0 |) is found to be smaller in H0B2C2 compared to DyB2C2. 
Moreover, resonant diffraction at the Ho L-edge has shown that the quadrupolar 
moment in H0B2C2 is smaller than in DyB2C2 [IS]- 

The temperature dependences recorded with SXRD, ND and XRD show 
noticeable differences due to their different structure factors (Fig. [6|. The structure 
factor of the (0 1 -y) reflection recorded with XRD is proportional to |24j 

^01 it oc [cos(2<?!)i) - cos(2(/i3)](Qoiri ) (3) 

where 0„ is the angle of the principal axis of the Ho orbital at site n with respect to the 
(1 0) axis in the ab plane and (Qqi^) is the expectation value of the Ho 4/ time-even 
multipolar moment which is proportional to a weighted sum of ^2 ^^'^ ^2 |24) . 
Thus the diffracted intensity observed with XRD is proportional to the orbital, or 
time-even, order of the Ho 4/ shell. In contrast, the intensity at (0 |) observed with 
ND is proportional to the magnetic dipole moment of the Ho 4/ shell, in particular 
to the magnetic moment components in the ab plane as those are perpendicular to 
the scattering vector. Figure [6] shows that the temperature dependence observed with 
these methods is distinct and this suggests that the temperature dependence of the 
orbital and magnetic order parameters are themselves distinct. Both order parameters 
shows a power law dependence / oc C'ip^ Y^'' with the magnetic critical exponent 
/3m ~ 0.35, while the orbital order increases more abruptly with Pq k, 0.2. The 
temperature dependence of the SXRD intensity recorded at the (0 0^) depends on 
energy, polarization and azimuthal angle and has the ND and XRD curves as lower 
and upper extremities, respectively (Figure [6]). This supports the claim that SXRD 
measures the interference of the magnetic and orbital order parameters of which the 
relative contribution depends on the energy, incident polarization and azimuthal angle. 

To illustrate this point. Figure [Tl] presents the energy dependence of the magnetic 
contributions, A\ and and orbital contribution, A\, to the scattered intensity 
if — 90° and tt incident radiation as well as f = 45° and a incident radiation. While 
in the first the magnetic part dominates for feature D, in the latter the orbital part 
dominates in feature B, consistent with the temperature dependencies presented in 
Fig. [6] when compared to neutron and x-ray diffraction respectively. 

Distinct temperature dependences within a single resonance profile has been 
reported in Lag.sSri 5Mn04 [5S]. At particular energies of the Mn L2^3-edges, the 
temperature dependence of the orbital order reflection (3; | 0) resembles that of the 
orbital order parameter, while at other energies the temperature dependence resembles 
that of the Jahn- Teller distortion. RXD on the (1 2) and (10|) reflections in 
TbB2C2, performed at the Tb L-edge, have shown two distinct critical exponents, 
tentatively assigned to the magnetic and orbital order parameters |27|. In this case, 
~ 0.3, and Pq « 0.45, suggesting that the orbital order is induced by (and possibly 
competes with) the magnetic order. Distinct timescales have also been observed in 
the paramagnetic and paraorbital state of DyB2C2, where it has been observed that 
the 4/ orbital moments fluctuate faster than the magnetic moments [26]. For H0B2C2 
it has been suggested [12] that the magnetic order induces the orbital order at Tg 
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Figure 9. (left) Fitted Parameters in 
the AFM + AFQ phases of H0B2C2. 
(right) Reconstruction of the angular 
dependence of the 4/ charge density 
described by the values of in the table. 



This value is set to 1 



due to the formation of a quasi doublet ground somewhere below T^v due to near- 
crossing of two crystal field levels in the presence of magnetic exchange. Orbital 
fluctuations can inhibit the AFQ order and inelastic neutron scattering performed 
on a H0B2C2 sample [18j noted a broadening of the peaks at Tq which has been 
attributed to strong orbital fluctuations at the critical temperature. In contrast to 
the case of TbB2C2, we find that those components of the scattered intensity whose 
contribution is largely orbital have the smaller critical exponent, and a consequently 
sharper temperature dependence. This indicates that the magnetic reorientation is 
at least partially induced by the orbital ordering. There are not enough data points 
from the neutron scattering data to attempt a proper measurement of the critical 
exponent, but an illustrative attempt has been made to convey the differences between 
the neutron and x-ray scattering data in Figure [TOj 

Remarkably, a finite intensity has been observed with SXRD up to T = 5 K 
while the intensity observed with both the XRD and ND is zero (Figure [6]). At 
present the reason for this is unclear. SXRD is sensitive to all directions of the 
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Temperature (K) 

Figure 10. temperature dependence of the neutron and x-ray intensities, 
normalized to 3.3 K, with illustrative power law fits / oc The power 

law curves correspond to /9q 0.2 and /3jv/ ~ 0.35 

magnetic moments, as well as to higher order time-odd multipole moments whereas 
ND measures the in plane dipole moments only. Thus a dipole moment along the c-axis 
or higher order time-odd multipole moments or can give rise to finite SXRD 
while ND and XRD remain zero. This implies a magnetic reorientation transition 
from the incommensurate long range magnetic order between Tjv and T = 5 K to 
commensurate magnetic order below T = 5 K to combined magnetic and quadrupolar 
order below Tg. However, the energy dependence observed with SXRD at T = 4.9 
K (not shown) displays the characteristics of combined orbital and magnetic order. 
Possibly, short range correlations associated with the transition at Tg may become 
pinned at the surface leading to an observable signal above Tq. SXRD probes the 
surface region of the sample (« 200 nm) while XRD and ND probe « 100 /im and 
the whole sample respectively. An additional possibility arises from the inherent 
uncertainty in measuring the sample temperature, and it cannot be ruled out that the 
intensity measured at 5 K with SRXD should coincide with the intensity measured at 
4.8 K with XRD and ND, although it is problematic to explain the presence of the 
jump in intensity at the lower temperature if this is the case. Therefore we cannot 
exclude that the SXRD signal observed above Tg is due to either a surface effect 
or the uncertainty in the absolute temperature. We emphasize that, if present, such 
effects cannot account for the differences observed in the temperature dependencies. 
Moreover, a difference in the measured ordering temperature of 0.2 K has a minor 
effect on the SXRD recorded at T = 3.4 K, which is well below Tg. 

6. Conclusion 

H0B2C2 exhibits an ordered phase (phase III) with both antiferromagnetic (AFM) and 
antiferro-orbital (AFQ) order below Tq = 4.8 K. Soft resonant x-ray Bragg diffraction 
experiments were performed at the Ho M4^5 edges and the energy dependence of 
the space group forbidden (00 1) reflection was measured in this low temperature 
AFM-I-AFQ phase. The energy dependence of the (00 i) reflection is modeled 
successfully by including an intra-atomic core-hole interaction parameterized by 
quadrupole and magnetic interactions. A pseudo-multiplet structure was included, 
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Energy (eV) Energy (eV) 

Figure 11. The energy dependence of the magnetic parts, A\ and Af, and 
orbital part, A|, of the (0 0^) reflection for (a) ip = 90° and n incident radiation 
and (b) ip = 45° and a incident radiation. B and D are labeled in accordance 
with Figure[4] Note that the orbital contribution is dominant for feature B in (b) 
while the magnetic contribution is dominant for feature D in (a), agreeing with 
the measured temperature dependences in Figure |6] 



allowing for the interference between oscillators that have different properties under 
time-reversal due to their having different excited-state lifetimes. The Ho 4/ multipole 
moments of rank 1 (dipole) to 6 (hexacontatetrapole) were determined, with the 
contributions from rank 1 and 3 time-odd multipoles found to be small compared 
to the dominant contribution of the hexadecapole (rank 4). 

Comparison of the temperature dependence of the orbital ordering reflections 
recorded with SXRD, XRD and ND shows that the magnetic and quadrupolar order 
parameters are distinct. The magnetic order exhibits a power-law dependence below 
Tq with a critical exponent that is smaller than that of the orbital order. We conclude 
the quadrupolar interaction is strong, but quadrupolar order only occurs when the 
magnetic order results in a near-crossing of the crystal field levels, producing a quasi 
doublet ground state. This is consistent with strong quadrupolar fluctuations present 
above Tg and the lock-in like transition of the orbitals at Tq. 
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Appendix A. The atomic scattering amplitude A^J,M and the reduced 
matrix elements R^{r,x) 

A^{J,M) is constructed from the structure factor of the chemical unit cell 4*^ and 
reduced matrix elements R^{r,x) according to 

<(J,M).(-l/-^^(2.+l)(_4 :)^"(^>-)*^(A.i) 

where the symbol in brackets is a 3jm symbol, with r — 0,1,..., 2 J, x = \K — 
r|, . . . , jif + r| and g + x and r + x are both even integers. The reduced matrix 
elements R^{r,x) are calculated for the Ho'^+ ground state of ^/g according to the 
relation found in the appendix of [H] ■ A Matlab program for calculating these reduced 
matrix elements for a valence configuration of n^^ or /-electron holes appropriate 
to either an El — El or E2 — E2 event is available from the corresponding authors 
upon request. 

Appendix B. Ho site symmetry and magnetic structure 

The structure factor is calculated using the magnetic structure presented in Figure [l] 
For simplicity the small c-axis canting described by |28| is excluded and we have the 



following structure factor of the atomic tensors. 

= {1 ~ e'^){{r;^) - (TZg)} (B.l) 

*f = -i (1.29(r|)' - 3.53(T|)") (B.2) 

= ~i (1.88i(rf )' - 1.32(Tf )") (B.3) 



where (. . .)' and (. . .)" represent respectively the real and imaginary parts of the 
atomic tensors. We note that the fitting process includes directly. 

Appendix C. Contributions to the scattering factor as a function of 
azimuthal angle 

The energy-independent component of the structure factor for contributions to the 
scattering is given (with the quantum numbers J and M of A^ omitted for the sake 



of brevity) as 

= -~2V5iAl cos(2(^) am{0) + \fhiX{ sm((p) cos{0) 

-\/3yl} cos(v5) cos(0) (C.l) 

fa,a = 2V5iAj sm{2ip) (C.2) 

/tt.ct = ~/(T,Tr (C.3) 

= VbiAl sin(2^)(l - cos(26l)) + VSA] am{ip) sm{2e) (C.4) 
which gives for Lp = 45° and 20 w 90° 

I, = ^\V5iAj~V3Al\'' + 20\Al\^ (C.5) 

I, = ^\V3Al - VbiAll^ + \V5iAl - Y^^ll' (C.6) 
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and for ip = 90° with 261 w 90° 



I, = l\2iAl+tAl\' 

I^ = 3\Al\^ + l\2iAl+iA'i\^ 



(C.7) 



(C.8) 
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